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There  is  increasing  interest  in  developing  solid  oxide  fuel  cells  (SOFC)  for  portable  applications.  For  these 
devices  it  would  be  convenient  to  directly  use  a  liquid  fuel  such  as  methanol  and  ethanol  rather  than  hydro¬ 
gen.  The  direct  utilization  of  alcohol  fuels  in  SOFC  involves  several  processes,  including  the  deposition  of 
carbon,  which  can  lead  to  irreversible  deactivation  of  the  fuel  cell.  Several  publications  have  addressed 
the  thermodynamic  analysis  of  the  reforming  of  methanol  (MeOH)  and  ethanol  (EtOH)  in  SOFC,  but  none 
have  considered  the  direct  utilization  of  these  fuels.  The  equilibrium  compositions,  the  carbon  deposition 
boundaries,  and  the  electromotive  forces  for  the  direct  utilization  and  partial  oxidation  of  methanol  and 
ethanol  in  SOFC  as  a  function  of  the  fuel  utilization  are  obtained  in  this  study.  In  addition,  the  minimum 
amounts  of  H20,  and  C02  for  direct  and  indirect  reforming  with  MeOH  and  EtOH  to  avoid  carbon  formation 
are  calculated. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  currently  considered  as  the  most  efficient  and 
clean  alternative  to  convert  the  chemical  energy  of  a  fuel  into 
electricity,  and  extensive  research  is  ongoing  to  improve  the  reli¬ 
ability  and  reduce  the  cost  of  these  electrochemical  devices  [1]. 
Most  fuel  cell  types  require  hydrogen  as  fuel  [1],  but  hydrogen  is 
expensive  and  difficult  to  store  if  compared  to  traditional  hydro¬ 
carbon  fuels.  Thus,  there  is  an  increasing  interest  in  developing 
fuel  cells  that  also  operate  with  traditional  fuels,  such  as  natural 
gas,  propane,  gasoline,  and  diesel.  Among  all  types,  solid  oxide  fuel 
cells  (SOFC)  offer  better  fuel  versatility;  in  fact,  synthesis  gas  is 
typically  used  to  fuel  SOFC  by  integrating  a  fuel-reforming  unit 
within  the  system  [2].  There  also  have  been  several  reports  on 
SOFC  operating  with  partial  oxidation  of  conventional  fuels  (e.g. 
propane  [3])  and  with  the  direct  utilization  of  natural  gas  and 
other  hydrocarbons  [4,5].  However,  the  utilization  of  synthesis  gas 
and,  in  particular,  the  direct  utilization  of  hydrocarbons  can  lead 
to  rapid  deactivation  by  sulfur  poisoning  and/or  carbon  formation 
on  the  traditional  SOFC  anode  [2]  (i.e.  a  composite  of  nickel  and 
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yttria-stabilized  zirconia,  Ni/YSZ).  In  the  case  of  direct  utilization, 
carbon  formation  can  be  particularly  severe,  causing  the  delamina¬ 
tion  of  the  Ni/YSZ  anode  [6]  that  rapidly  and  irreversibly  destroys 
the  cell. 

Although  SOFC  are  intended  mainly  for  medium-large  station¬ 
ary  applications,  where  a  reforming  unit  can  be  easily  integrated, 
there  is  an  increasing  interest  in  developing  SOFC  for  portable 
applications,  such  as  portable  power  generation  units  and  auxiliary 
power  units  in  military  and  transportation  uses  [7].  For  these  appli¬ 
cations,  it  would  be  ideal  to  use  a  liquid  fuel,  which  is  easy  to  store 
and  transport,  and  to  operate  the  cell  in  the  direct  utilization  or 
partial  oxidation  modes,  which  allows  a  reduction  of  the  size,  com¬ 
plexity,  and  cost  of  the  system  [2].  Methanol  (MeOH)  and  ethanol 
(EtOH)  have  recently  received  much  attention  as  liquid  fuels  partic¬ 
ularly  as  alternative  ‘energy- vectors’  for  the  future  [8].  Both  of  these 
fuels  have  a  very  low  content  of  impurities  (<0.2%),  are  oxygenated, 
which  drastically  reduces  the  risks  of  sulfur-poisoning  and  carbon 
formation,  and  are  candidate  fuels  for  portable  SOFC  applications. 

The  direct  utilization  of  alcohols  in  SOFC  involves  many  reac¬ 
tion  processes  in  both  the  gas  phase  and  at  the  anode,  including: 
(1)  pyrolysis  in  the  gas-channels,  (2)  catalytic  decomposition  on 
the  anode  surface,  (3)  electrochemical  oxidation  in  proximity  of 
the  triple-phase  boundary,  (4)  steam  and  dry  reforming  on  the 
anode  surface  and/or  oxidation  of  the  fuel,  (5)  water-gas  shift  reac¬ 
tion  on  the  anode  surface,  and  (6)  carbon  deposition  and  removal. 
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Nomenclature 

CDB 

carbon  deposition  boundary 

E° 

equilibrium  potential 

EMF 

electromotive  force 

ES;_j 

equilibrium  selectivity  for  species  i  and  j 

EtOH 

ethanol 

Eyt 

equilibrium  yield  for  species  i 

F 

Faraday  constant  (96485.34  C  mol-1 ) 

MeOH 

methanol 

ft  fuel-in 

molar  flow  rate  of  alcohol  fuel  at  the  inlet 

ftfuel-eq 

molar  flow  rate  of  alcohol  fuel  at  the  outlet  at  equi¬ 
librium 

ft-i-eq 

molar  flow  rate  of  species  i  at  the  outlet  at  equilib¬ 
rium 

ft02~echem 

molar  flow  rate  of  oxygen  due  to  electrochemical 
oxidation 

ft(D  2-leak 

molar  flow  rate  of  oxygen  due  to  air  leaks 

ft02-P0X 

molar  flow  rate  of  oxygen  due  to  POX 

OCP 

open  circuit  potential 

POX 

partial  oxidation 

RFh2o 

steam-reforming  factor 

RFco2 

dry-reforming  factor 

SOFC 

solid  oxide  fuel  cell 

Uf 

fuel  utilization 

Uf-app 

apparent  fuel  utilization 

Xmax 

maximum  fuel  conversion 

Z 

number  of  electrons  produced  per  alcohol  molecule 

The  extent  of  each  reaction  is  related  to  the  flow  conditions,  tem¬ 
perature  gradients,  and  overall  current.  Carbon  deposition,  which 
can  deactivate  and/or  mechanically  destroy  the  anode,  will  happen 
whenever  there  is  a  positive  driving  force  for  the  deposition  reac¬ 
tions.  According  to  the  affinity  concept  (see  [9]),  this  driving  force 
is  in  the  first  approximation  proportional  to  the  difference  between 
the  actual  and  the  equilibrium  activities  for  the  reacting  species. 

Typically  the  SOFC  operating  conditions  (e.g.  steam  to  carbon 
ratio)  are  chosen  to  completely  avoid  carbon  formation.  Although 
SOFC  may  not  operate  at  equilibrium,  thermodynamic  analysis  of 
gas-phase  compositions  and  carbon  formation  provides  a  useful 
guideline.  Several  studies  on  the  equilibrium  for  dry  and  steam 
reforming  have  been  reported  [10-14],  and  few  studies  on  the  inter¬ 
nal  reforming  and  partial  oxidation  of  methanol  and  ethanol  [15,16]. 
The  focus  of  these  publications  was  on  the  identification  of  the  min¬ 
imum  FI20,  C02  and  air  requirements  to  prevent  carbon  deposition 
for  external  or  internal  reforming  of  several  hydrocarbon  fuels,  but 
no  attention  was  given  to  the  variations  of  the  equilibrium  con¬ 
ditions  with  the  degree  of  oxidation  of  the  fuel.  Only  one  study 
addressed  the  variation  of  the  carbon  deposition  boundaries  as  a 
function  of  the  current  density  for  the  direct  utilization  of  methane 
in  SOFC  [17].  Cairns  et  al.  [18,19]  provided  a  general  treatment  of  the 
equilibrium  for  low  temperature  (<500  K)  hydrocarbon  fuel  cells.  In 
fact,  understanding  the  variations  of  equilibrium  composition  and 
electromotive  force  (EMF)  as  a  function  of  the  fuel  utilization  is 
fundamentally  important  for  an  SOFC  operating  in  the  direct  uti¬ 
lization  and  partial  oxidation  modes.  Thus,  the  work  of  Cairns  et  al. 
[18,19]  was  used  as  a  guide  for  the  study  herein. 

The  objective  of  the  work  reported  in  this  paper  is  to  obtain 
equilibrium  compositions,  carbon  deposition  boundaries,  and  elec¬ 
tromotive  forces  for  the  direct  utilization  and  partial  oxidation  of 
MeOFI  and  EtOH  in  SOFC.  To  reach  this  objective  the  pyrolysis, 
steam-reforming,  and  dry-reforming  conditions  were  studied  with 
and  without  the  presence  of  a  current.  The  minimum  amounts 


of  H20,  C02  and  02  to  avoid  carbon  deposition  for  SOFC  fuelled 
directly  or  indirectly  with  MeOH  and  EtOH  are  also  reported.  All  val¬ 
ues  were  obtained  using  a  nonstoichiometric  formulation  approach 
and  calculations  were  performed  with  the  software  Chemical  Equi¬ 
librium  with  Applications  (CEA)  developed  by  NASA. 

2.  Theory  and  calculations 

Equilibrium  compositions  at  the  anode  of  an  SOFC  were  calcu¬ 
lated  by  assuming  perfectly  mixed  isothermal  reactor  conditions 
for  the  anodic  compartment  (i.e.  homogeneous  temperature  and 
pressure)  with  fuel-gas  fed  at  a  constant  rate.  The  amounts  of  each 
species  present  at  equilibrium  in  the  gaseous  and  solid  phases  were 
obtained  in  the  temperature  range  between  773  and  1173  K,  at  inter¬ 
vals  of  25  K,  and  for  a  total  pressure  of  1  bar.  The  type  of  carbon 
that  forms  at  an  anode  is  complex  and  depends  on  many  factors 
[20].  For  the  majority  of  the  present  calculations,  graphite  was  con¬ 
sidered  as  the  only  allotropic  species  for  solid  carbon  present  at 
equilibrium.  Some  calculations  were  performed  to  assess  the  impli¬ 
cations  of  including  amorphous  carbon,  rather  than  graphite,  in  the 
equilibrium  calculations. 

At  constant  temperature  and  pressure  the  Gibbs  free  energy  for 
a  system  at  equilibrium  is  at  a  global  minimum.  The  basic  criterion 
for  the  establishment  of  chemical  equilibrium  is  given  by: 


(dCSYs)T,p  =  0  (1) 

This  criterion  was  used  for  the  determination  of  the  equilibrium 
compositions  in  the  anodic  compartment.  There  are  two  different 
approaches  to  solve  this  problem.  The  first,  called  stoichiometric 
formulation  [21],  essentially  requires  the  solution  of  a  system  of 
nonlinear  equations  obtained  using  the  stoichiometric  equations 
of  the  chemical  reactions  considered  to  obtain  the  minimum  of 
Eq.  (1).  The  second  approach,  called  nonstoichiometric  formula¬ 
tion,  minimizes  the  Gibbs  free  energy  as  a  function  of  the  moles 
of  a  large  collection  of  species  using,  for  example,  the  method  of 
Lagrange  multipliers,  and  does  not  require  the  assumption  of  any 
chemical  reactions.  The  solutions  obtained  using  the  two  proce¬ 
dures  are  equivalent,  provided  that  a  complete  set  of  reactions 
is  used  in  the  first  approach  [21].  The  equilibrium  compositions 
presented  in  this  study  were  calculated  numerically  using  the  soft¬ 
ware  CEA  (Chemical  Equilibrium  with  Applications)  developed  by 
NASA.1  CEA  implements  the  nonstoichiometric  formulation,  and 
uses  a  database  of  thermodynamic  properties  containing  over  2000 
species,  including  solids,  liquids  and  gases.  The  results  obtained 
using  CEA  for  a  few  basic  cases  and  the  implications  of  consid¬ 
ering  amorphous  carbon  rather  than  graphite  were  verified  using 
the  commercial  software  for  thermodynamic  calculations  FactSage 
5.4, 2  which  also  utilizes  a  nonstoichiometric  method.  In  all  cases 
checked,  the  equilibrium  results  obtained  by  CEA  and  FactSage  were 
essentially  identical. 

The  equilibrium  compositions  plotted  in  this  work  are  expressed 
as  maximal  fractional  conversion,  equilibrium  yield  and  selectivity 
according  to  the  following  definitions  [22]: 


^max  — 


Eyt  = 


ftfuel  -in  ftfuel- eq 
ftfuel-in 


ni-( 


i-eq 


ftfuel-in  Kfuel-eq 


ESh  = 


li-eq 

nj-eq 


(2) 

(3) 

(4) 


1  www.grc.nasa.gov/WWW/CEAWeb/. 

2  www.factsage.com. 
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where  Xmax  is  the  fuel  conversion  at  equilibrium,  Eyt  is  the  equi¬ 
librium  yield  for  species  i,  ES ,-_j  is  the  equilibrium  selectivity  for 
species  i  relative  to  species  j,  rijuej_in  and  rifuei_eq  are  the  initial  and 
equilibrium  molar  flow  rates,  respectively,  of  the  fuel  and  xi[_eq  is 
the  equilibrium  molar  flow  rate  for  species  i. 

When  an  SOFC  is  producing  a  current  there  is  a  net  input  of 
oxygen  ions  (O2-)  to  the  anodic  compartment  that  changes  the 
equilibrium  conditions  and  so  calculations  were  performed  for 
various  current  loads.  In  the  case  of  direct  utilization,  the  new 
equilibrium  composition  can  be  predicted  considering  that  the  oxy¬ 
gen  inflow  due  to  electrochemical  oxidation  is  related  to  the  fuel 
utilization  by  Eq.  (5): 

fto2-echem  =  0-25  Z  Uf  Tlfuei-in  (5) 

where  z  is  the  number  of  electrons  produced  by  the  oxidation  of 
1  alcohol  molecule  (i.e.  6  electrons  for  MeOH  and  12  electrons  for 
EtOH),  Uq2  -echem  is  the  molar  flow  rate  of  oxygen,  and  Uf  is  the  fuel 
utilization.  The  fuel  utilization  (Uf)  is  related  to  the  overall  current 
produced  by  the  cell  (/tot)  by: 


/tot  =  zFUfrifuei_in  (6) 

where  F  is  the  Faraday  constant.  In  the  case  of  partial  oxidation, 
there  is  an  additional  input  of  02.  Similarly,  the  presence  of  air  leaks 
further  increases  the  02  inflow,  while  the  presence  of  N2  dilutes  the 
feed.  The  total  oxygen  input  to  the  anodic  compartment  is  given  by 
the  sum  of  these  inflows,  and  the  apparent  fuel  utilization  is  given 
by: 


Uf-app  ~  Uf  + 


no2-POX  +  n o2-leak 

0.25  zrifuei_in 


(7) 


where  n02_pox  and  n02_ieak  are  the  moles  of  oxygen  admitted  with 
the  fuel  (in  the  case  of  POX)  and  through  any  leaks,  respectively. 
In  this  study,  the  equilibrium  compositions  were  calculated  as  a 
function  of  the  fuel  utilization  for  mixtures  of  alcohol-fuel  and  02 
according  the  ratio  given  by  Eq.  (5). 

The  carbon  deposition  boundaries  (CDB)  were  obtained  by  cal¬ 
culating  the  equilibrium  composition  over  a  range  of  alcohol-fuel 
to  H20,  C02,  or  02  ratios.  The  carbon  formation  region  is  identified 
by  the  conditions  for  which  the  amount  of  graphite  present  at  equi¬ 
librium  is  larger  than  1 E-6  mole  fractions.  Likewise,  the  carbon  free 
region  corresponds  to  conditions  for  which  graphite  at  equilibrium 
is  smaller  than  IE-6  mole  fractions. 

The  electromotive  force  was  calculated  using  the  Nernst 
equation  including  all  species  present  at  equilibrium  (including 
graphite)  at  the  anode,  and  only  air  at  the  cathode  (p02  =  0.206PTOt)- 
The  heat  capacities,  enthalpies  and  entropies  of  formation  used 
were  those  tabulated  by  the  database  DIPPR  [23].  The  activity  for 
graphite  was  assumed  equal  to  unity,  while  the  fugacity  of  the 
gaseous  species  was  assumed  equivalent  to  the  equilibrium  partial 
pressures. 

In  a  real  cell,  the  maximum  fuel  utilization  that  can  be  reached 
for  a  given  fuel  flow  rate  depends  on  the  actual  cell  performance 
(i.e.  i-V  curve)  according  to  Eq.  (8): 


Uf-  max 


^lim 

zFllfuel-m 


(8) 


where  /iim  is  the  limiting  or  the  maximum  current  for  the  cell.  The 
maximum  current  achievable  with  the  SOFC  depends  on  several 
factors,  e.g.  ohmic  and  mass  transfer  limitations,  electrochemical 
overpotential  of  anode  and  cathode.  Eqs.  (7)  and  (8)  can  be  used 
to  verify  whether  an  experimental  SOFC  is  operating  outside  the 
carbon  formation  region  during  experiments  of  direct  utilization 
with  alcohols. 


3.  Results  and  discussion 

3.1.  Equilibrium  compositions  with  methanol  feed 

The  equilibrium  composition  for  the  pyrolysis  of  MeOH  as  a 
function  of  temperature  is  shown  in  Fig.  1 .  In  all  figures,  the  symbols 
represent  points  at  which  calculations  were  performed,  while  the 
lines  represent  the  trend.  The  only  species  present  at  equilibrium 
are  hydrogen  (H2),  carbon  monoxide  (CO),  methane  (CH4),  carbon 
dioxide  (C02),  water  (H20)  and  graphite  (C).  The  equilibrium  con¬ 
version  of  MeOH  is  100%  at  all  temperatures  (i.e.  MeOH  is  entirely 
decomposed).  The  presence  of  C  is  predicted  over  the  entire  tem¬ 
perature  range.  The  equilibrium  concentration  of  C  decreases  with 
temperature,  from  16.6  mol%  at  773  K  to  only  0.2mol%  at  1173  K. 
Within  the  operational  conditions  typical  of  SOFC,  i.e.  973-1073  K, 
the  concentration  of  C  changes  from  7.5  to  1.9  mol%,  respectively, 
indicating  that  the  ‘driving  force’  for  the  carbon  deposition  reaction 
decreases  with  temperature.  In  practice,  carbon  deposition  may  be 
highly  dependent  on  kinetics. 

As  shown  in  Fig.  lb,  the  gas-phase  is  composed  mostly  by  H2 
(approximately  62  vol%)  and  CO  (approximately  25  vol%).  CH4,  C02 
and  H20  become  less  prevalent  at  higher  temperatures.  The  H2 
to  CO  ratio  (H2/CO)  approaches  2  with  increasing  temperature.  At 
1073 1<  the  H2/CO  is  2.1,  indicating  that  the  direct  decomposition 
to  H2  and  CO  (i.e.  CH3OH  =  2H2  +  CO)  is  the  predominant  reaction 
(i.e.  most  negative  Gibbs  free  energy  change  and  largest  equilib¬ 
rium  constant)  among  all  the  possible  reactions  involving  MeOH. 
This  result  can  be  verified  by  applying  the  stoichiometric  method 
for  the  decomposition  reaction.  As  a  result,  the  gas-phase  corn- 


800  900  1000  1100 

T  emperature  /K 


T  emperature  /  K 

Fig.  1.  Equilibrium  composition  (a)  and  gas-phase  concentrations  (b)  as  a  function 
of  temperature  for  the  reaction  of  MeOH  pyrolysis.  The  species  at  equilibrium  are 
(•)  CH4,  (O)  CO,  (▼)  C02,  (□)  H2,  (■)  H20,  and  (□)  C  (graphite). 
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Table  1 

Minimum  reforming  factors  required  to  avoid  carbon  formation  in  the  steam-  and  dry- reforming  of  methanol,  and  the  resulting  H2  yield,  CH4  yield  and  H2/CO  ratios. 


T( K)  MeOH  steam-reforming  MeOH  dry-reforming 


RFh2o 

Eyn2 

Eye  h4 

ESh2/co 

RFco2 

Eyu2 

Eye  h4 

ESh2/co 

773 

0.848 

0.827 

0.526 

12.186 

38.000 

0.503 

0.013 

0.206 

823 

0.828 

1.153 

0.422 

7.180 

14.500 

0.720 

0.028 

0.332 

873 

0.704 

1.446 

0.309 

4.559 

5.500 

0.973 

0.050 

0.533 

923 

0.499 

1.650 

0.209 

3.197 

1.900 

1.252 

0.078 

0.872 

973 

0.288 

1.762 

0.134 

2.514 

0.614 

1.511 

0.087 

1.322 

1023 

0.140 

1.827 

0.086 

2.200 

0.205 

1.700 

0.072 

1.681 

1073 

0.060 

1.873 

0.056 

2.070 

0.072 

1.817 

0.052 

1.865 

1123 

0.023 

1.908 

0.037 

2.020 

0.025 

1.887 

0.036 

1.946 

1173 

0.007 

1.934 

0.025 

2.002 

0.007 

1.927 

0.025 

1.981 

position  in  this  case,  can  be  well  estimated  by  only  considering 
the  MeOH  decomposition  reaction.  Species  such  as  formaldehyde 
(HCHO)  and  ethylene  (C2H4)  are  not  found  at  equilibrium  because 
they  are  unstable  relative  to  CH4,  CO,  and  C02.  In  fact,  HCHO  is  an 
important  intermediate  of  decomposition  and  is  found  when  the 
pyrolysis  of  MeOH  is  not  complete  [24].  The  results  shown  in  Fig.  1 
are  representative  of  the  equilibrium  state  in  the  anode  compart¬ 
ment  of  a  SOFC  at  OCP  (i.e.  zero  current)  when  the  fuel  consists  of 
pure  methanol. 

Equilibrium  calculations  were  performed  to  identify  the  con¬ 
ditions  at  which  carbon  deposition  can  be  avoided  for  steam  and 
dry  reforming  of  MeOH;  the  results  are  summarized  in  Table  1.  The 
steam-reforming  factor  (RFH2o)  and  dry-reforming  factor  (RFCo2) 
are  defined  as  the  molar  ratios  of  H20  to  MeOH  and  C02  to  MeOH, 
respectively,  required  to  prevent  carbon  (as  graphite)  deposition. 
The  equilibrium  compositions  for  the  given  reforming  factors  are 
reported  in  terms  of  equilibrium  yields  (EyH2 ,  Eyc h4  )  and  selectiv¬ 
ity  (ESh2/co ,  ESh2/ch4  ) as  defined  in  Eqs.  (3 )  and  (4).  Note  that  the 
equilibrium  conversion  for  both  steam  and  dry  reforming  is  100% 
over  the  entire  temperature  range. 

For  the  steam-reforming  of  MeOH  the  minimum  value  of  the 
reforming  factor  that  should  be  used  to  avoid  carbon  formation  is 
0.060  at  1073 1<  (i.e.  6  mol  of  H20  per  100  mol  of  MeOH).  Consider¬ 
ing  that  dry  MeOH  is  hygroscopic,  it  is  evident  that  the  use  of  highly 
pure  feedstock  is  paramount  if  the  direct  utilization  of  this  fuel  is  to 
be  studied.  In  fact,  if  the  liquid  MeOH  feedstock  used  to  fuel  a  SOFC 
at  1073  K  contains  more  than  2.6  vol%  of  H20,  the  mode  of  opera¬ 
tion  is  direct  internal  reforming  rather  than  direct  utilization.  The 
reforming  factor  decreases  with  temperature,  in  accordance  with 
the  results  previously  shown  for  MeOH  pyrolysis  (Fig.  1 ).  Consider¬ 
ing  the  equilibrium  compositions  corresponding  to  the  minimum 
reforming  factors,  for  temperatures  higher  than  1073  K,  ESH2/co  is 
approximately  2,  indicating  that  MeOH  has  formed  almost  exclu¬ 
sively  H2  and  CO. 

At  temperatures  lower  than  923  K,  the  dry-reforming  of  MeOH 
requires  high  ratios  of  C02  to  MeOH  (between  1.9  and  38)  to  avoid 
carbon  formation.  At  773  K,  an  equimolar  mixture  of  MeOH  and  C02 
produces  20  mol%  more  carbon  at  equilibrium  than  MeOH  pyrol¬ 
ysis.  Conversely,  at  temperatures  greater  than  1073  K  both  FyH2 
and  ESh2/co  approach  2,  indicating  that  MeOH  has  essentially  only 
formed  H2  and  CO,  similar  to  the  result  for  steam-reforming  of 
MeOH. 

The  equilibrium  compositions  were  calculated  as  a  function  of 
the  fuel  utilization  considering  an  inlet  stream  composed  of  dry 
MeOH  with  oxygen  from  the  cathode  in  the  proportions  given  by  Eq. 
(5).  The  region  of  carbon  deposition  shown  in  Fig.  2a  is  identified  by 
the  minimum  fuel  utilization  for  which  graphite  is  absent  (<lE-6 
mole  fraction)  from  the  equilibrium  mixture  at  a  given  temperature. 
The  CDB  shifts  towards  lower  fuel  utilization  as  the  temperature 
increases.  The  CDB  as  a  function  of  the  fuel  utilization  shown  in 
Fig.  2a  also  can  be  used  to  determine  the  minimum  amount  of 


air  required  to  avoid  coke  formation  for  the  partial  oxidation  of 
methanol  (MeOH-POX).  In  fact,  by  using  Eq.  (7)  it  is  possible  to  cal¬ 
culate  the  necessary  flow  rate  of  oxygen  to  bring  the  system  outside 
the  carbon  formation  region.  For  example,  a  SOFC  fuelled  directly 
with  MeOH  at  973  K  is  outside  the  carbon  formation  region  if  the 
fuel  utilization  is  13.25%.  For  POX  operation,  this  condition  must 
be  achieved  also  at  OCP  (i.e.  7  =  0  A;  Uf=  0%).  Therefore,  it  is  neces¬ 
sary  to  artificially  increase  the  fuel  utilization  (i.e.  Uf_app  =  13.25%) 
by  introducing  air  with  the  fuel.  Assuming  that  air  leaks  are  negligi¬ 
ble  (i.e.  Uq2 -leaks  =  0  mol  s_1 ),  from  Eq.  (7)  the  amount  of  02  to  add 
with  MeOH  is  no2-pox  =  Uf-app  x  0-25  x  z  x  nfuei  =  0.198  x  n^e/,  or 
nair  =  0.964rifuei,  mols-1.  The  minimum  02  to  MeOH  ratios  for 
MeOH-POX  as  a  function  of  the  temperature  are  summarized  in 
Table  2.  The  corresponding  equilibrium  compositions  and  poten¬ 
tials,  as  a  function  of  the  temperature,  using  the  minimum  amount 
of  air  required  to  avoid  coking  are  shown  in  Fig.  2b. 
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Fig.  2.  Carbon  formation  region  as  a  function  of  the  temperature  (a)  for  a  SOFC 
fuelled  with  dry  MeOH.  Equilibrium  composition  potential  (£°)  for  the  partial  oxi¬ 
dation  of  MeOH  (b)  at  the  minimum  air  to  fuel  ratio  to  avoid  coking.  The  species  at 
equilibrium  are  (•)  CH4,  (O)  CO,  (▼)  C02,  (A)  H2,  and  (■)  H20. 
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Table  2 

Minimum  oxygen  to  fuel  ratios  (no2/n^e/)  for  MeOH  and  EtOH  partial  oxidation  to 
avoid  carbon  formation. 


T(K) 

( n02  lnfuei ) 

MeOH 

EtOH 

773 

1.1E+00 

2.3E+00 

823 

9.2E-01 

2.0E+00 

873 

6.8E-01 

1.7E+00 

923 

4.1E-01 

1.3E+00 

973 

2.0E-01 

9.5E-01 

1023 

8.4E-02 

7.2E-01 

1073 

3.3E-02 

6.0E-01 

1123 

1.2E-02 

5.4E-01 

1173 

3.4E-03 

5.2E-01 

Finally,  the  anode  compositions  and  electromotive  forces  for 
the  direct  utilization  of  dry  MeOH  as  a  function  of  fuel  utiliza¬ 
tion  at  three  temperatures  (1073,  973  and  873  K)  were  calculated 
(Fig.  3).  These  three  graphs  have  some  common  features.  As  the  fuel 
utilization  increases  the  concentrations  of  H20  and  C02  increase, 
while  the  concentration  of  H2  decreases.  The  concentration  of  CO 
increases  with  increasing  fuel  utilization  within  the  carbon  forma¬ 
tion  region,  while  it  decreases  when  the  fuel  utilization  is  higher 
than  the  CDB.  The  H2/CO  does  not  change  significantly  with  the 
fuel  utilization.  At  1073  K,  the  H2/CO  is  approximately  2  even  up  to 
80%  fuel  utilization.  The  EMFs  at  0%  fuel  utilization  (i.e.  the  equi¬ 
librium  potential,  £°)  are  1.084, 1.057,  and  1.041  V  at  1073,  973  and 
873  K,  respectively.  As  expected,  the  EMF  decreases  with  increas¬ 
ing  fuel  utilization  and  the  rate  of  this  decrease  is  higher  at  higher 
temperatures.  For  example,  at  1073  K  the  EMF  decrease  going  from 
0%  to  80%  fuel  utilization  is  0.181  V,  while  at  873 1<  the  decrease 
is  0.088  V. 

3.2.  Equilibrium  compositions  with  ethanol  feed 

The  equilibrium  composition  for  the  pyrolysis  of  EtOH  as  a 
function  of  temperature  is  shown  in  Fig.  4.  As  for  MeOH,  the 
EtOH  conversion  at  equilibrium  is  always  100%,  and  the  only  sta¬ 
ble  species  are  H2,  CO,  CH4,  C02,  H20,  and  C.  Species  such  as 
acetaldehyde  (CH3CHO),  acetic  acid  (CH3COOH),  ethylene  (C2H4) 
and  formaldehyde  (HCHO)  are  not  predicted,  but  are  likely  to  be 
important  intermediates  in  the  actual  process.  Carbon  deposition 
is  predicted  over  the  entire  temperature  range  considered.  The 
amount  of  carbon  at  equilibrium  decreases  with  the  temperature 
stabilizing  to  approximately  20  mol%  above  1123  K,  at  which  EtOH  is 
entirely  decomposed  to  H2,  CO  and  C.  The  ratio  of  these  species  is  in 
accordance  with  the  reaction  C2H5OH  =  3H2  +  CO  +  C,  such  that  only 
this  reaction  is  required  to  determine  the  equilibrium  composition 
above  1123  K. 

Table  3  contains  the  reforming  factors,  yields,  and  selectivities 
for  ethanol  steam-  and  dry-reforming  at  temperatures  between 
773  and  1173  K.  For  the  steam-reforming  of  EtOH,  the  minimum 
value  of  the  reforming  factor  that  should  be  used  to  avoid  car¬ 
bon  formation  is  1.125  at  1073  K  (i.e.  112.5  mol  of  H20  per  100  mol 
of  EtOH).  As  shown  earlier,  the  minimum  steam-reforming  fac¬ 
tor  for  MeOH  is  only  0.060  (20  times  lower),  which  makes  direct 
utilization  almost  indistinguishable  from  direct  internal  reforming 
in  the  case  of  methanol.  A  steam-reforming  factor  equal  to  1.125 
corresponds  to  an  EtOH-H20  liquid  mixture  containing  25.7  vol% 
H20,  which  is  more  than  5  times  larger  than  in  the  azeotropic 
mixture  and  more  than  3  times  than  in  the  typical  commercial  alco¬ 
hol  formulation.  The  reforming  factor  corresponding  to  the  CDB 
decreases  with  temperature,  and  approaches  1  for  temperatures 
above  1123  K  in  accordance  with  the  stoichiometry  of  EtOH  reform  - 


Fig.  3.  Equilibrium  composition  and  electromotive  force  for  a  SOFC  operated  with 
dry  MeOH  at  (a)  1073  K,  (b)  973  K,  and  (c)  873 1<.  The  species  at  equilibrium  are  (•) 
CH4,  (O)  CO,  (▼)  C02,  (A)  H2,  (■)  H20,  and  (n)  C  (graphite). 


Table  3 

Minimum  reforming  factors  required  to  avoid  carbon  formation  in  the  steam-  and 
dry-  reforming  of  ethanol,  and  the  resulting  H2  yield,  CH4  yield  and  H2/CO  ratios. 


T  (K)  EtOH  steam-reforming  EtOH  dry- reforming 


RFh2o 

Eyn2 

Eye  h4 

ESh2/co 

RFco2 

£yH2 

Eye  h4 

ESh2/co 

773 

2.709 

1.660 

1.051 

12.232 

81 

0.756 

0.014 

0.146 

823 

2.670 

2.314 

0.841 

7.209 

33 

1.075 

0.028 

0.224 

873 

2.416 

2.898 

0.617 

4.568 

14 

1.449 

0.050 

0.333 

923 

2.051 

3.341 

0.407 

3.244 

6.25 

1.846 

0.072 

0.478 

973 

1.583 

3.530 

0.267 

2.519 

3.2 

2.210 

0.074 

0.633 

1023 

1.269 

3.700 

0.176 

2.149 

1.85 

2.516 

0.069 

0.785 

1073 

1.125 

3.787 

0.100 

2.090 

1.35 

2.712 

0.053 

0.882 

1123 

1.046 

3.818 

0.073 

2.020 

1.15 

2.780 

0.031 

0.925 

1173 

1.014 

3.869 

0.050 

2.002 

1.07 

2.893 

0.027 

0.964 
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Fig.  4.  Equilibrium  composition  (a)  and  gas-phase  composition  (b)  as  a  function  of 
the  temperature  for  the  reaction  of  EtOH  pyrolysis.  The  species  at  equilibrium  are  (•) 
CH4,  (O)  CO,  (▼)  C02,  (A)  H2>  (■)  H20,  and  (□)  C  (graphite).  Equilibrium  composition 
potential  (£°)  for  the  partial  oxidation  of  EtOH  (b)  at  the  minimum  air  to  fuel  ratio 
to  avoid  coking.  The  species  at  equilibrium  are  (•)  CH4,  (O)  CO,  (▼)  C02,  (A)  H2,  and 
(■)H20. 

ing  (i.e.  CH3CH20H  +  H20  =  4H2 +2CO).  Consistently,  the  yield  of 
hydrogen  (£yH2)  increases  with  temperature  and  approaches  a 
value  of  4,  while  the  hydrogen  to  carbon  monoxide  ratio  (ESH2/co) 
decreases  with  temperature  and  approaches  a  value  of  2.  The  quan¬ 
tity  of  CH4  present  at  equilibrium  decreases  with  temperature, 
consistent  with  the  conversion  to  H2  and  CO  being  the  dominant 
reaction. 

As  for  MeOH,  the  dry-reforming  of  EtOH  requires  relatively 
high  reforming  factors  to  avoid  carbon  formation  at  temper¬ 
atures  lower  than  973  K.  Also  in  this  case,  the  presence  of 
C02  promotes  carbon  formation  at  773 1<  for  reforming  factors 
smaller  than  2  if  compared  to  EtOH  pyrolysis.  The  reforming 
ratio  decreases  with  temperature  and  approaches  1  for  temper¬ 
atures  higher  than  1123  K  as  in  the  stoichiometry  of  the  EtOH 
dry-reforming  reaction  (i.e.  CH3CH20H  +  C02  =  3H2  +3CO).  Sim¬ 
ilarly,  EyH 2  increases  with  temperature  approaching  3,  while 
ESh2/co  approaches  1.  The  equilibrium  yield  for  CH4  is  always 
low  relative  to  that  for  H2  and  CO  and  reaches  a  maximum  at 
approximately  923  K. 

Equilibrium  compositions  were  calculated  as  a  function  of  fuel 
utilization  considering  an  inlet  stream  composed  of  dry  EtOH  and 
oxygen  transported  through  the  electrolyte  in  the  proportions  given 
by  Eq.  (5).  From  the  results  obtained,  the  region  of  carbon  deposi¬ 
tion  was  identified  as  shown  in  Fig.  5a.  The  required  fuel  utilization 
required  to  avoid  carbon  formation  decreased  with  temperature. 
For  example,  at  973  K  the  minimum  fuel  utilization  to  avoid  coking 
is  28.7%,  while  at  1073  K  it  is  19.5%.  Unlike  the  case  with  MeOH, 


T  emperature  IK 

Fig.  5.  Carbon  formation  region  as  a  function  of  temperature  for  a  SOFC  fuelled  with 
dry  EtOH. 


the  CDB  line  does  not  approach  zero  at  high  temperatures  (com¬ 
pare  with  Fig.  2).  The  equilibrium  composition  and  potential,  as  a 
function  of  temperature,  in  the  case  of  partial  oxidation  of  ethanol 
(EtOH-POX)  using  the  minimum  amount  of  air  required  to  avoid 
coking  is  shown  in  Fig.  5b.  The  minimum  02  to  fuel  ratio  for  the 
POX  of  EtOH  as  a  function  of  the  temperature  are  summarized  in 
Table  2.  Because  the  extent  of  carbon  formation  is  greater  with 
EtOH  than  MeOH,  the  minimum  levels  of  02  are  between  approxi¬ 
mately  2  (at  773  K)  and  approximately  150  (at  1173  K)  times  larger 
for  EtOH. 

The  equilibrium  compositions  and  EMF  for  the  direct  utilization 
of  EtOH  as  a  function  of  fuel  utilization  at  1073,  973  and  873 1<  are 
shown  in  Fig.  6a-c,  respectively.  The  trends  of  the  stable  species 
in  Fig.  6  are  similar  to  those  seen  for  MeOH  in  Fig.  3.  Specifically, 
with  increasing  fuel  utilization  the  concentrations  of  H20  and  C02 
increase,  while  the  concentration  of  H2  decreases.  The  concentra¬ 
tion  of  CO  initially  increases,  reaches  a  maximum  at  the  CDB,  and 
then  decreases  with  fuel  utilization.  The  H2/CO  ratio  decreases  for 
fuel  utilizations  greater  than  the  CDB.  The  concentration  of  CH4 
decreases  significantly  with  temperature.  At  873  K,  the  maximum 
CH4  concentration  is  14.6  vol%,  while  at  1073  K  the  maximum  con¬ 
centration  is  2.4  vol%  (also  see  Fig.  4a).  In  both  cases,  these  maxima 
occur  at  0%  fuel  utilization.  The  EMFs  at  0%  fuel  utilization  are 
1.103,  1.068,  and  1.052  V  at  1073,  973  and  873  K,  respectively.  The 
EMF  decreases  with  increasing  fuel  utilization,  decreasing  approxi¬ 
mately  0.225  V  when  going  from  0%  to  80%  fuel  utilization  at  1073  K. 
As  for  MeOH,  the  change  in  EMF  with  fuel  utilization  decreases  with 
decreasing  temperatures. 
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Fig.  6.  Equilibrium  composition  and  electromotive  force  for  a  SOFC  operated  with 
dry  EtOH  at  (a)  1073  K,  (b)  973  K,  and  (c)  873  K.  The  species  at  equilibrium  are  (•) 
CH4,  (O)  CO,  (▼)  C02,  (A)  H2,  (■)  H20,  and  (□)  C  (graphite). 


3.3.  Effect  of  using  thermodynamic  data  for  amorphous  carbon 

Several  studies  have  shown  that  the  type  of  carbon  deposited 
in  SOFC  depends  on  the  type  of  hydrocarbon  used  as  fuel  [2],  on 
the  anode  material  [25],  and  possibly  also  on  other  factors,  such 
as  aging  and  history  of  the  sample  [20].  Thermodynamic  measure¬ 
ments  indicate  that  graphite  is  more  stable  than  amorphous  carbon. 
Unfortunately,  there  are  only  a  few  studies  reporting  the  enthalpy 
and  entropy  of  formation  for  amorphous  carbon  [26-28],  and  there 
is  no  agreement  on  the  enthalpy  and  entropy  difference  between 
amorphous  carbon  relative  to  graphite  [29,30].  According  to  [28] 
the  enthalpy  of  formation  of  amorphous  carbon  is  15.400  kj  mol-1 
higher  than  that  of  graphite,  while  the  entropy  is  4.1  JI<-1  mol-1 
higher,  compatible  to  the  higher  disorder  of  the  amorphous  struc¬ 
ture.  The  standard  Gibbs  free  energy  for  the  transformation  of 
amorphous  carbon  into  graphite  was  found  to  be  -11.054  kj  mol-1 . 
In  another  study  [26],  the  enthalpy  difference  varied  significantly 
depending  on  the  degree  of  amorphousness  of  the  carbon,  but  the 
standard  Gibbs  free  energy  for  the  transformation  of  amorphous 


carbon  into  graphite  was  in  all  cases  smaller  than  -3.220  kj  mol-1 . 
The  effect  on  the  equilibrium  composition  of  using  thermodynamic 
data  for  amorphous  carbon  rather  than  graphite  was  investigated 
for  the  Boudard  reaction  (2CO  =  C(s) +  C02)  and  for  the  methane 
decomposition  (CH4  =  C(S)  +  2H2 )  using  the  thermodynamic  data  for 
amorphous  carbon  proposed  in  [28]. 

The  equilibrium  constant  for  the  Boudard  reaction  calculated 
using  the  standard  thermodynamic  data  at  773  K  for  graphite  is 
6.7  times  higher  than  that  obtained  using  the  data  for  amorphous 
carbon,  while  at  1173  K  the  equilibrium  constant  is  2.9  times  higher 
for  graphite.  The  ratios  of  the  equilibrium  constant  for  the  methane 
decomposition  reaction  vary  in  a  similar  fashion.  Thus,  the  effect 
of  using  the  values  for  amorphous  carbon  instead  of  the  values  for 
graphite,  is  a  decrease  of  the  moles  of  carbon  (C(s))  present  at  equi¬ 
librium  for  both  MeOH  and  EtOH.  If  both  graphite  and  amorphous 
carbon  are  included  in  the  calculation  of  the  equilibrium  composi¬ 
tion,  the  moles  of  amorphous  carbon  are  negligible  in  all  cases,  and 
the  resulting  composition  is  identical  to  that  obtained  using  only 
the  data  for  graphite.  This  result  was  expected  since  the  energy 
content  of  amorphous  carbon  is  higher  that  graphite.  In  addition,  it 
is  expected  that  amorphous  carbon  should  transform  into  graphite 
over  time  (i.e.  aging),  as  this  reaction  is  exothermic. 

4.  Conclusions 

According  to  the  results  of  the  equilibrium  calculations  pre¬ 
sented  in  this  paper,  H2  and  CO  are  the  main  products  of  MeOH 
pyrolysis.  The  presence  of  C  is  predicted  over  the  temperature  range 
considered  but  the  equilibrium  concentration  decreases  signifi¬ 
cantly  with  temperature,  becoming  less  than  2  mol%  above  1073  K. 
At  these  temperatures,  the  equilibrium  constant  for  the  reaction 
CH3OH  =  2H2  +  CO  is  higher  than  any  other  decomposition  pathway 
such  that  only  this  reaction  needs  to  be  considered  to  estimate  the 
equilibrium  composition.  The  presence  of  small  quantities  of  H20 
in  the  (liquid)  MeOH  used  to  fuel  the  SOFC  (e.g.  2.6  vol%  for  a  SOFC 
operating  at  1073  K)  should  be  sufficient  to  suppress  carbon  forma¬ 
tion  at  equilibrium;  therefore,  the  purity  of  the  feed  is  important  in 
experimentation  to  ensure  that  the  SOFC  is  operated  in  the  direct 
utilization  mode.  Carbon  formation  also  becomes  thermodynami¬ 
cally  unfavourable  when  the  SOFC  is  producing  current.  The  min¬ 
imum  values  of  fuel  utilizations  for  which  graphite  is  not  present 
at  equilibrium  were  identified.  For  MeOH,  the  fuel  utilization  value 
at  the  carbon  formation  boundary  is  2.1%  at  1073  K,  and  higher  fuel 
utilizations  are  needed  to  operate  in  the  carbon  free  region. 

For  EtOH  pyrolysis,  the  main  products  at  equilibrium  are  H2,  CO 
and  C.  The  reaction  C2H5OH  =  3H2 +  CO  +  C  dominates  the  equilib¬ 
rium  of  decomposition  above  1123  K.  For  EtOH  steam-reforming, 
relatively  large  amounts  of  H20  must  be  added  to  pure  EtOH  in 
order  to  prevent  carbon  formation  (e.g.  25.7  vol%  for  a  SOFC  oper¬ 
ating  at  1023  K).  The  carbon  deposition  boundary  was  calculated  as 
a  function  of  the  fuel  utilization.  At  1073  K,  the  minimum  fuel  uti¬ 
lization  for  a  SOFC  fuelled  with  dehydrated  EtOH  to  operate  outside 
the  carbon  formation  region  is  19.5%. 

The  use  of  thermodynamic  data  for  amorphous  carbon  rather 
than  graphite  in  the  calculation  of  the  thermodynamic  equilib¬ 
rium  results  in  a  decrease  in  the  amount  of  solid  carbon  present 
at  equilibrium.  Because  the  energy  content  of  amorphous  carbon 
is  significantly  higher  than  graphite,  the  use  of  both  amorphous 
carbon  and  graphite  in  the  calculations  results  in  the  same  final 
equilibrium  composition  as  using  only  graphite. 

In  an  operating  SOFC  there  might  be  significant  deviations  from 
equilibrium.  One  such  deviation,  for  instance,  is  related  to  the 
incomplete  conversion  of  the  alcohols.  However,  the  equilibrium 
conditions  described  in  this  paper  are  indicative  of  the  final  state 
towards  which  the  system  may  evolve. 
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